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Abstract — An adaptive randomized distributed space-time cod- 
ing (DSTC) sclieme and algoritlims are proposed for two-liop 
cooperative MIMO networks. Linear minimum mean square 
error (MMSE) receivers and an amplify-and-forward (AF) coop- 
eration strategy are considered. In the proposed DSTC scheme, a 
randomized matrix obtained by a feedbaclt channel is employed 
to transform the space-time coded matrix at the relay node. 
Linear MMSE expressions are devised to compute the parameters 
of the adaptive randomized matrix and the linear receive filter. 
A stochastic gradient algorithm is also developed to compute 
the parameters of the adaptive randomized matrix with reduced 
computational complexity. We also derive the upper bound of 
the error probability of a cooperative MIMO system employing 
the randomized space-time coding scheme first. The simulation 
results show that the proposed algorithms obtain significant 
performance gains as compared to existing DSTC schemes. 

I. Introduction 

Multiple-input and multiple-output (MIMO) communication 
systems employ multiple collocated antennas at both the 
source node and the destination node in order to obtain the 
diversity gain and combat multi-path fading in wireless links. 
The different methods of space-time coding (STC) schemes, 
which can provide a higher diversity gain and coding gain 
compared to an uncoded system, are also utilized in MIMO 
wireless systems for different numbers of antennas at the 
transmitter and different conditions of the channel. Coop- 
erative MIMO systems, which employ multiple relay nodes 
with antennas between the source node and the destination 
node as a distributed antenna array, apply distributed diversity 
gain and provide copies of the transmitted signals to improve 
the reliability of wireless communication systems [l]. Among 
the links between the relay nodes and the destination node, 
cooperation strategies, such as Amplify-and-Forward (AF), 
Decode-and-Forward (DF), and Compress-and-Forward (CF) 
lis and various distributed STC (DSTC) schemes in 13], (U 
and |l20l| can be employed. 

The utilization of a distiibuted STC (DSTC) at the relay 
node in a cooperative network, providing more copies of the 
desired symbols at the destination node, can offer the system 
diversity gains and coding gains to combat the interference. 
The recent focus on the DSTC technique lies in the delay- 
toleration code design and the full-diversity schemes design 
with the minimum outage probability. In (|5], the distributed 
delay-tolerant version of the Golden code ^ is proposed, 
which can provide full-diversity gain with a full coding rate. 
An opportunistic DSTC scheme with the minimum outage 
probability is designed for a DF cooperative network and 
compared with the fixed DSTC schemes in l?). An adaptive 
distributed-Alamouti (D-Alamouti) STBC design is proposed 



in fE\ for the non-regenerative dual-hop wireless system which 
achieves the minimum outage probability. DSTC schemes for 
the AF protocol are discussed in |l9]-|[l0l. In M, the GABBA 
STC scheme is extended to a distributed MIMO network with 
full-diversity and full-rate, while an optimal algorithm for 
design of the DSTC scheme to achieve the optimal diversity 
and multiplexing tradeoff is derived in 110|. 

The performance of cooperative networks using different 
strategies has been widely discussed in the literature. In fTTl, 
an exact pairwise error probability of the D-Alamouti STBC 
scheme is derived according to the position of the relay node. 
In 1 12J . a bit error rate (BER) analysis of the distributed- 
Alamouti STBC scheme is proposed. The difference between 
these two works lies in the different cooperative schemes 
considered. A maximum likelihood (ML) detection algorithm 
for a MIMO relay system with DF protocol is derived in 
lT3l with its performance analysis as well. The symbol error 
rate and diversity order upper bound for the scalar fixed- 
gain AF cooperative protocol are given in [14|. The use of 
single-antenna relay nodes and the DF cooperative protocol 
is the main difference in scenario between lfT4l and this 
work. An STC encoding process is implemented at the source 
node in |18|, which decreases the output of the system and 
increases the computational complexity of the decoding at 
the destination node. In IJ9J . the BER upper bound is given 
without a STC scheme at the relay node. 

In this paper, we propose an adaptive randomized distributed 
space-time coding scheme and algorithms for a two-hop coop- 
erative MIMO relaying system with the AF protocol and linear 
MMSE receivers. We focus on how the randomized matrix 
affects the DSTC during the encoding and how to optimize 
the parameters in the matrix. It is shown that the use of a 
randomized matrix benefits the performance of the system by 
lowering the upper bound compared to using traditional STC 
schemes. Linear MMSE expressions are devised to compute 
the parameters of the adaptive randomized matrix and the 
linear receive filter. Then an adaptive optimization algorithm 
is derived based on the MSE criterion, with the stochastic 
gradient (SG) algorithm in order to reduce the computational 
complexity of the optimization process. The updated random- 
ized matrix is transmitted to the relay node through a feedback 
channel that is assumed in this work error free and delay free . 
The upper bound pairwise error probability of the randomize d- 
STC schemes (RSTC) in a cooperative MIMO system which 
employs multi-antenna relay nodes with the AF protocol is 
also analyzed. 

The paper is organized as follows. Section II introduces 
a two-hop cooperative MIMO system with multiple relays 
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applying the AF strategy and the randomized DSTC scheme. 
In Section III the proposed MMSE expressions and the SG 
algorithm for the randomized matrix are derived, and the 
analysis of the upper bound of pairwise error probability using 
the randomized D-STC is shown in Section IV. Section V 
focus on the results of the simulations and Section VI leads 
to the conclusion. 

II. Cooperative MIMO System Model 

The communication system under consideration is a two- 
hop cooperative MIMO system employing multiple relay 
nodes and communicating over channels from the source node 
to the relay nodes and the destination node, and from the 
relay nodes to the destination nodes as shown in Fig. 1. A 
modulation scheme is used in our system to generate the 
transmitted symbol vector at the source node and all 
nodes have iV antennas. There are relay nodes the employ 
the AF cooperative strategy as well as a DSTC scheme. The 
system broadcasts symbols from the source to relay nodes 
as well as to the destination node in the first phase. The 
symbols are amplified and re-encoded at each relay node 
prior to transmission to the destination node in the second 
phase. We consider only one user at the source node in our 
system that has N Spatial Multiplexing (SM)-organized data 
symbols contained in each packet. The received symbols at 
the k ~ th relay node and the destination node are denoted 
as rsB^k ^^'^ fSD, respectively, where k = l,2,...,nr. The 
received symbols rsB.^ arc amplified before mapped into an 
STC matrix. We assume that the synchronization at each node 
is perfect. The received symbols at the destination node and 
each relay node are described by 



rsR^ m 
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i = 1,2, 



Fk\i]s[i\+nsR,[i\ 
- H[i\s[{\ + nsoii], 
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(1) 
(2) 



where the x 1 vector nsR^.[i] and n5£)[i] denote the zero 
mean complex circular symmetric additive white Gaussian 
noise (AWGN) vector generated at each relay and the destina- 
tion node with variance a^. The transmitted symbol vector s[i] 
contains N parameters, s[i] — [si[i], S2[i], SNii]], which 
has a covariance matrix ii^[s[z]s-'^[j]] = a^J, where E[-] 
stands for expected value, (•) denotes the Hermitian operator, 
(Tg is the signal power which we assume to be equal to 1 
and / is the identity matrix. Fk[i] and H[i] are the N x N 
channel gain matrices between the source node and the kth 



relay node, and between the source node and the destination 
node, respectively. 

After processing and amplifying the received vector 
f'SB.ki'^] at the kth relay node, the signal vector ssi?fcH — 
^i^fcuWl-F + ■'T'SflfcW) can be obtained and then for- 
warded to the destination node. The amplified symbols in 
^SRk W will be re-encoded by a A^xTDSTC scheme M(s[i]) 
and then multiplied by an x randomized matrix in 
ll24ll . then forwarded to the destination node. The relationship 
between the k — th relay and the destination node can be 
described as 



Rr,d[^ = Gfc[z]CH[z]Mfl,z5W + NR,D[i], 
— 1 , 2 , . . . , , 



(3) 



where the A x T matrix Mn^^oli] is the DSTC matrix 
employed at the relay nodes whose elements are the amplified 
symbols in S57?j,[i]. The A x T received symbol matrix 
-RflfccM ill © can be written as an AT x 1 vector rji^ali] 
given by 

rR^D [i] = 0\eq^ \i] Geq^ [i] SSR^ \i] + D W , (4) 

where the block diagonal AT x AT matrix 0\eqk [*] denotes 
the equivalent randomized matrix and the AT x A matrix 
Geq^ [i] Stands for the equivalent channel matrix which is the 
DSTC scheme Af (§[?]) combined with the channel matrix 
Gi?j.£)[i]. The AT x 1 equivalent noise vector nflj.£i[i] gen- 
erated at the destination node contains the noise parameters 
in MB.^nii]- After rewriting R^.^Dii] we can consider the 
received symbol vector at the destination node as a N{nr + 1) 
vector with two parts, one is from the source node and another 
one is the superposition of the received vectors from each 
relay node. Therefore, the received symbol vector for the 
cooperative MIMO system can be written as 



r I 



= DD[i]sD[i\ + noli], 

(5) 

where the (T + 1)A x (n,. + 1)A^ block diagonal matrix 
denotes the channel gain matrix of all the links in 
the network which contains the N x N channel coefficients 
matrix H[i] between the source node and the destination 
node, the A^T x A equivalent channel matrix Geq^ [i] for 
k = 1, 2, ...,nr between each relay node and the destination 
node. The (n^ + 1)A^ x 1 noise vector n£)[i] contains the 
received noise vector at the destination node and the amplified 
noise vectors from each relay node, which can be modeled 
as additive white Gaussian noise (AWGN) with zero mean 
and covariance matrix (7^{l+ \\ 9ieqk['i]Geq^^[^]AIi^D[^] 11%)!, 

where || X \\f= \J^^{X" ■ X) = 
Frobenius norm. 
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III. Design of Linear MMSE Receivers and 
Randomized Matrices 

In this section, we design an adaptive linear MMSE receive 
filter and an MMSE randomized matrix for use with the 
proposed DSTC scheme. An adaptive SG algorithm |21| for 
determining the parameters of the randomized matrix with 
reduced complexity is also devised. The DSTC scheme used at 



the relay node employs an MMSE randomized matrix, which 
is computed at the destination node and obtained by a feedback 
channel and processes the data symbols prior to transmission 
to the destination node. 

A. Optimization Method Based on the MSE Criterion 

Let us consider the MMSE design of the receive filter and 
the randomized matrix according to the optimization problem 



[T^H,91egM] = arg mm E 



w 



Hi 



where r[i] is the received symbol vector at the destination node 
which contains the randomized matrix to be optimized. If we 
only consider the received symbols from the relay node, the 
received symbol vector at the destination node can be derived 

as 

r[i] = DD[i]sD[i] + noii] 

= meq[l]Geq[l]A[i]F[i]s[i]+^eq[i]Geq[i]A[i]nSR[i] 

+ riRD \i] 

(6) 

where C[i] is an NT x N matrix that contains all the 
complex channel gains and the amplified matrix assigned 
to the received vectors at the relay node, and the noise 
vector is a Gaussian noise with zero mean and variance 
(7^(1+ II 9legHGeg[i]A[i] \\p). We can then recast the 
optimization as 



arg min E 



(7) 

By expanding the righthand side of (|7]i and taking the gradient 
with respect to W* [i] and equating the terms to zero, we can 
obtain the linear MMSE receive filter 

W[i] = {E [r[^r"[^])''E[r[i]s"[i]\ , (8) 

where the first term denotes the inverse of the auto-correlation 
matrix and the second one is the cross-correlation matrix. 
Define r — C[i]s[i] + C[i]nsFi^, then the randomized matrix 
can be calculated by taking the gradient with respect to OV [i] 
and equating the terms to zero, resulting in 

= (w"[i]iE [r[i]r"[i]])W[{\y' E[s[i]r"[i]] W[i], 

(9) 

where E [r[i]r^[i]\ is the auto-correlation of the space- 
time coded received symbol vector at the relay node, and 
E [s[i\r"[i\] is the cross-correlation. The expression above 
requires a matrix inversion with a high computational com- 
plexity. 

B. Adaptive Randomized Matrix Optimization Algorithm 

In order to reduce the computational complexity and achieve 
the optimal performance, an adaptive randomized matrix op- 
timization (ARMO) algorithm based on an SG algorithm 
is devised. The MMSE problem is derived in (|7]i, and the 
MMSE filter matrix can be calculated by ^ first during the 
optimization process. The simple ARMO algorithm can be 



obtained by taking the instantaneous gradient term of (|7]l with 
respect to the randomized matrix lHeq*[i], which is given by 



VE 



= -is{i] - W"[i]r[i])s"[i]C"\i]W[i] 
= -e[i]s"[i]C"[i]W[i], 

(10) 

where e[i] stands for the detected error vector After we 
computing ( fTOl ). the ARMO algorithm can be obtained by 
introducing a step size into an SG algorithm to update the 
result until the convergence is reached as given by 

m[i + i] = m[i] + fi{e[i]s"[i]c"[i]w[i]), (11) 

where p stands for the step size of the ARMO algorithm. The 
complexity of calculating the randomized matrix is 0{2N), 
which is much less than that of the calculation method derived 
in As mentioned in Section I, the randomized matrix will 
be sent back to the relay nodes via a feedback channel which 
is assumed to be error- free in this work. However, in practical 
circumstances, the errors caused by the broadcasting and the 
diversification of the feedback channel with time changes will 
affect the accuracy of the received randomized matrix at the 
relay nodes. 

IV. Probability of Error Analysis 

In this section, the upper bound of the pairwise error 
probability of the system employing the randomized DSTC 
will be derived. As we mentioned in the first section, the 
randomized matrix will be considered in the derivation as 
it affects the performance by reducing the upper bound of 
the pairwise error probability. For the sake of simplicity, we 
consider a 2 by 2 MIMO system with 1 relay node, and the 
direct link is ignored in order to concentrate on the effect of 
the randomized matrix. The expression of the upper bound is 
also stable for the increase of the system size and the number 
of relay nodes. 

Consider an x STC scheme we use at the relay 
node with L codewords. The codeword is transmitted 
and decoded to another codeword at the destination node, 
where i = 1,2,...,L. According to |25|, the probability of 
error can be upper bounded by the sum of all the probabilities 
of incorrect decoding, which is given by 



(12) 



Assuming the codeword C is decoded at the destination 
node and we know the channel information perfectly at the 
destination node, we can derive the pairwise error probability 

as 

P(C^ ^ I 91) 
= P(|| R} - GOIC^ III - II i?^ - GOIC^ |||> I Oleq) 

-P(|| -OleqGeqFs^ \\% 

-\\r^ -9iegGeqFs^ \\l>Q\Oleq), 

(13) 

where F and Geq stand for the channel coefficient matrix 
between the source node and the relay node, and between 



the relay node and the destination node, respectively. The 
randomized matrix is denoted by 9ieq- Define H — G^qF, 
which stands for the total channel coefficients matrix. After 
the calculation, we can transfer the pairwise error probability 
expression in ( fTST l to 

P(Cl ^ I y\,q) = P(|| ^eqH{s^ - ||2,< Y), (14) 

where r = Tr{'n}" 'y\eqH{s^-s^) + {^^qH{s^-s^))"'n}), 
and denotes the noise vector at the destination node with 
zero mean and covariance matrix cr^(l+ || £HegGeg 
By making use of the Q function, we can derive the error 
probability function as 



P(Cl ^ I ^eq) = Q II "neqHis^ 



-s') II 



where 



-= / exp -— du, 



(15) 
(16) 



and 7 is the received SNR at the destination node assuming 
the transmit power is equal to 1 . 

In order to obtain the upper bound of P(C^ ^ \ *He,) 
we expand the formula II yieqH{s^ — s'^) |||-. Let AgU he 
the eigenvalue decomposition of {s^ — s^)^(s^ — s^), where 
U is a Hermitian matrix and contains all the eigenvalues 
of the difference between two different codewords and 
s^. Let Ay{V stand for the eigenvalue decomposition of 
{0leqHU)"9leqHU, whcrc V is a random Hermitian matrix 
and A(K is the ordered diagonal eigenvalue matrix. Therefore, 
the probability of error can be written as 



Oleq) = 



NT N 



n.rn | 



m—l n—1 



(17) 

where ^n,m is the {n, m) — th element in V, and A<h,^ and 
As„ are eigenvalues in At^ and A^, respectively. According 
to ||251 . a good upper bound assumption of the Q function is 
given by 

Cl{x)<]^e^. (18) 

Thus, we can derive the upper bound of pairwise error prob- 
ability for a randomized STC scheme as 



P(Ci^C2|Ole,)<7Texp 



NT N \ 

^ E E ^9^™-^S,ll^":"l|^ 1 

m—l n — 1 / 

(19) 

while the upper bound of the error probability expression for 
a traditional STC is given by 

/ NT N \ 

With comparison of (fT9] l and ( l20l i. it is obvious to note that 
the eigenvalue of the randomized matrix is the difference, 
which suggests that employing a randomized matrix for a STC 
scheme at the relay node can provide an improvement in BER 
performance. 



V. Simulations 

The simulation results are provided in this section to assess 
the proposed scheme and algorithms. The cooperative MIMO 
system considered employs an AF protocol with the Alamouti 
STBC scheme |25| using QPSK modulation in a quasi-static 
block fading channel with AWGN, as described in Section 11. 
The bit error ratio (BER) performance of the ARMO algorithm 
is assessed. The simulation system with 1 relay node and each 
transmitting and receiving node employs 2 antennas. In the 
simulations, we define both the symbol power at the source 
node and the noise variance for each link as equal to 1. 

The upper bounds of the D- Alamouti and the randomized D- 
Alamouti derived in the previous section are shown in Fig. 2. 
The theoretical pairwise error probabilities provide the largest 
decoding errors of the two different coding schemes and as 
shown in the figure, by employing a randomized matrix at the 
relay node it decreases the decoding error upper bound. The 
comparison of the simulation results in BER performance of 
the R-Alamouti and the D-Alamouti indicates the advantage 
of using the randomized matrix. 

The proposed ARMO algorithm is compared with the 
SM scheme and the traditional RSTC algorithm using the 
distributed- Alamouti (D- Alamouti) STBC scheme in |20| with 
rir = 1 relay nodes in Fig. 3. The number of antennas 
= 2 at each node and the effect of the direct Unk are 
considered. The results illustrate that without the direct link, 
by making use of the STC or the RSTC technique, a significant 
performance improvement can be achieved compared to the 
spatial multiplexing system. The RSTC algorithm outperforms 
the STC-AF system, while the ARMO algorithm can improve 
the performance by about 3dB as compared to the RSTC 
algorithm. With the consideration of the direct link, the results 
indicate that the cooperative diversity order can be increased, 
and using the ARMO algorithm achieves an improved per- 
formance with 2dB of gain as compared to employing the 
RSTC algorithm and 3dB of gain as compared to employing 
the traditional STC-AF algorithm. 

The simulation results shown in Fig. 4 illustrate the con- 
vergence property of the ARMO algorithm. The SM, D- 
Alamouti and the randomized D-Alamouti algorithms obtain 
nearly flat performance in BER as the utilization of fixed 
STC scheme and the randomized matrix. The SM scheme 
has the worst performance due to the lack of coding gains, 
while the D-Alamouti scheme can provide a significant perfor- 
mance improvement in terms of the BER improvement, and 
by employing the randomized matrix at the relay node the 
BER performance can decrease further when the transmission 
circumstances are the same as that of the D-Alamouti. The 
ARMO algorithm shows its advantage in a fast convergence 
and a lower BER achievement. At the beginning of the 
optimization process with a small number of samples, the 
ARMO algorithm achieves the BER level of the D-Alamouti 
one, but with the increase of the received symbols, the ARMO 
algorithm achieves a better BER performance. 

VI. Conclusion 

We have proposed an adaptive randomized matrix optimiza- 
tion (ARMO) algorithm for the randomized DSTC using a 
linear MMSE receive filter at the destination node. The pair- 
wise error probability of introducing the randomized DSTC in 



a cooperative MIMO network with the AF protocol has been 
derived. The simulation results illustrate the advantage of the 
proposed ARMO algorithm by comparing it with the coop- 
erative network employing the traditional DSTC scheme and 
the fixed randomized STC scheme. The proposed algorithm 
can be used with different distributed STC schemes using the 
AF strategy and can also be extended to the DF cooperation 
protocol. 
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